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TOWARDS USING OPTIC AL/NIR PHOTOMETRY TO MEASURE THE 

TEMPERATURE OF O STARS 
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RESUMEN 

Tradicionalmente se ha dicho que no es posible usar fotometria visible/IR para medir la temperatura de las 
estrellas O. En esta contribucion describimos los pasos necesarios para veneer los obstaculos existentes hasta 
la fecha y presentamos nucstros rcsultados prcliminarcs para cl caso dc cxtincion baja. 

ABSTRACT 

It has been traditionally stated that it is not possible to use optical/NIR photometry to measure the tem- 
peratures of O stars. In this contribution we describe the steps required to overcome the hurdles that have 
prevented this from happening in the past and we present our preliminary results for the low-extinction case. 
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1. BACKGROUND 

Hummer et al. (1988) declared two decades ago 
that "for hot stars (Tcff > 30 000 K), methods based 
on the integrated continuum flux are completely un- 
reliable discriminators of the effective temperature" . 
This statement has been subsequently reformulated 
to express that optical/NIR broad- or medium-band 
photometry cannot be used to determine the intrin- 
sic properties of O stars. However, the data used by 
Hummer et al. (1988), as the authors themselves rec- 
ognize, had observational uncertainties of 0.02 mag- 
nitudes in E{B — V) and their flux calibration was 
only accurate to within 10% in the UV and 3% in 
the optical/NIR. Therefore, it is a legitimate ques- 
tion to ask whether optical/NIR can indeed be used 
to measure 0-star temperatures if the precision and 
accuracy of our data are better than those amounts. 
In this contribution we present the steps required to 
reduce the uncertainties to a 1% level. 

2. THE INTRINSIC COLORS OF O STARS 

2.1. Filter properties and zero-point determination 

The first required step is to calibrate the pho- 
tometry by: (a) accurately measuring the reference 
SED (e.g. Vega), (b) determining the zero point that 
fine-tunes the correspondence between the observed 
and the reference SEDs, and (c) correctly defining 
the total (atmosphere -|- telescope -I- filter -I- detec- 
tor) sensitivity curve of the system, (a) and (b) 
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are needed in order to avoid systematic errors be- 
tween different filters while (c) is needed to avoid 
"color terms" within a given filter. Some recent pa- 
pers have managed to produce significant improve- 
ments in all of the above. Bohlin (2006) has used 
a combination of STIS spectrophotometry and Ku- 
rucz models to obtain a new SED for Vega that elim- 
inates previous discrepancies between the results of 
different groups. Several other papers (Cohen et al. 
2003, Mai'z Apellaniz 2005, 2006a, 2006b, Holberg 
& Bergeron 2006) have analyzed different filter sys- 
tems (Johnson, Stromgren, 2MASS. . . ) to test the 
validity of the published sensitivity curves and to 
calculate the zero points. Interestingly, independent 
results for the same filter agree to within 1%, with 
the only outstanding exception of 2MASS Kg (Mafz 
Apellaniz 2006b). 

2.2. Atmosphere models 

The second required step is to be able to use at- 
mosphere models as inputs for the calculation of the 
intrinsic colors of O stars. The last five years have 
seen a significant increase in the detail of such mod- 
els. The two most significant additions have been 
the inclusion of line blanketing by heavy elements 
and wind effects, the former being the one that has 
the largest effect in the optical and the latter the 
dominant one in the K band and longer wavelengths 
(when mass-loss rates arc high, see Martins & Plez 
2006). There are now grids for early B (Lanz & 
Hubeny 2006) and O (Lanz & Hubeny 2003, Mar- 
tins et al. 2005) stars that cover the Tgg — \ogg — Z 
space needed to derive synthetic broad-band colors 
for most hot stars. Furthermore, the synthetic colors 
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Fig. 1. Difference between photometric and spectro- 
scopic temperature determinations for six low-extinction 
O stars using different SED families as inputs. 

from different models agree to within ~ 1%, lending 
credence to their reliability and making them appro- 
priate to be used as inputs for temperature determi- 
nations of early-type stars. 

3. TEMPERATURE DETERMINATION 
TECHNIQUES: CHORIZOS 

The traditional method for calculating Toff from 
photometry has been through color-color diagrams 
and Q-parameter determinations. That strategy has 
several problems: [a] the slope of an extinction tra- 
jectory in a color-color plot depends on the initial 
SED; [b] the slope is not constant, either, as a func- 
tion of the amount of extinction; and [c] it only al- 
lows the simultaneous use of information from two 
colors, hence being especially sensitive to system- 
atic errors in the calibration and to deviations from 
the expected unextincted SEDs. All of the above 
can easily introduce considerable systematic errors in 
the determination of stellar temperatures but those 
problems arc minimized when one uses a Bayesian 
code like CHORIZOS (Maiz Apellaniz 2004), which 
can process information from many filters simulta- 
neously without using the constant-slope Q approx- 
imation. This constitutes the third required step to- 
wards determining Tes for O stars from photometry. 

Taking into account all of the previous consider- 
ations and first test of the use of CHORIZOS 
to determine 0-star temperatures, we are working 
on an analysis of optical/NIR photometry of low- 



extinction {E{B -V) < 0.15) O stars with CHORI- 
ZOS. Our preliminary results for six late-type O 
stars using Stromgren and 2MASS photometry fit- 
ting temperature, extinction, and extinction law si- 
multaneously are shown in Fig. 1. Using either 
TLUSTY or CMFGEN models as input SEDs yields 
photometric temperatures with random uncertain- 
ties of 1000-2000 K and systematic differences of 
< 1 000 K with respect to the spectroscopic temper- 
atures (based on the Martins et al. 2005 calibration) . 
Kurucz SEDs yield similar random uncertainties but 
a larger systematic difference. 

4. EXTINCTION LAW 
The above results for low extinction arc of lit- 
tle practical value for Galactic stars because most 
O stars in the Milky Way are moderately or heav- 
ily extinguished, with an extinction law that varies 
for different sightlines. Accurately accounting for ex- 
tinction is the last required step in the process and, 
arguably, the hardest one to fulfill. The most com- 
monly used family of optical/NIR Galactic extinc- 
tion laws is that of Cardelli et al. (1989). How- 
ever, that paper was based on an analysis of a small 
(29) sample of stars with only moderate extinction 
and used an unphysical seventh degree polynomial in 
1/A for the optical region. Hence, it is not possible 
to dcrcdden observed photometry with the required 
accuracy for our goals (< 1%) using the Cardelli et 
al. (1989) extinction law for E{B - V) > 0.5. For 
that reason we believe that the time has come for 
the calculation of an alternative sightline-dependent 
Galactic optical/NIR extinction law. Some work has 
been recently done in the NIR (Moore et al. 2005, 
Nishiyama et al. 2006) but the optical has not been 
revised yet. We plan to do so in the incoming year 
using the data in the Galactic O star catalog as in- 
put. 
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